ENUMERATING THE BASINS OF INFINITY OF CUBIC 

POLYNOMIALS 

LAURA DE MARCO AND AARON SCHIFF 

Abstract. We study the dynamics of cubic polynomials restricted to their basins 
of infinity, and we enumerate topological conjugacy classes with given combinatorics. 



1. Introduction 

Let / : C ^ C be a polynomial of degree 3 with complex coefficients. Its basin of 
infinity is the open invariant subset 

X{f) = {zeC: r{z) ^ oo as n ^ oo}. 

In this article, we examine combinatorial topological-conjugacy invariants of the re- 
stricted dynamical system / : X{f) — > X{f) and count the number of possibilities 
for each invariant. We apply results of |DP1] which show that we can use these in- 
variants to classify topological conjugacy classes of pairs (/, X(/)) within the space 
of cubic polynomials. Moreover, when / is in the shift: locus, meaning that both of its 
critical points lie in X{f), these invariants classify conjugacy classes of polynomials 
/ : C ^ C. 

Specifically, we implement an algorithm which counts topological conjugacy classes 
of cubic polynomials of generic level N, defined by the condition that 

here {ci, €2} is the set of critical points of the polynomial / and 

Gfiz)= hm ;^log+|r(z)| 

is the escape-rate function. These generic cubic polynomials are precisely the struc- 
turally stable maps in the shift locus jMcH] (see also jDP2] l We examine the growth 
of the number of these stable conjugacy classes as — > cx). 

We begin with an enumeration of the Branner-Hubbard tableaux (or equivalently. 



the Yoccoz r-functions) of length A^, as introduced in |BH2j : see Theorem 2.2 Using 



the combinatorics of tableaux, we provide an algorithm for computing the number of 



truncated spines (introduced in |DP1] ) for each r-function; see Theorem 3.1 Finally, 



we apply the procedure of |DP1] to count the number of generic topological conjugacy 
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Level 


Tau sequences 


Trees 


Truncated spines 


Conjugacy classes 


Run time 


1 


1 


1 


1 


1 


0.000 


2 


2 


2 


2 


2 


0.078 


3 


4 


4 


4 


4 


0.062 


4 


8 


8 


8 


8 


0.063 


5 


16 


18 


18 


19 


0.093 


6 


33 


42 


42 


46 


0.079 


7 


69 


103 


105 


118 


0.078 


8 


144 


260 


270 


318 


0.093 


9 


303 


670 


718 


881 


0.094 


10 


641 


1753 


1939 


2480 


0.125 


11 


1361 


4644 


5312 


7084 


0.156 


12 


2895 


12433 


14719 


20374 


0.266 


13 


6174 


33581 


41161 


59061 


0.547 


14 


13188 


91399 


115856 


172016 


1.141 


15 


28229 


250452 


328098 


503018 


2.453 


16 


60515 


690429 


933719 


1475478 


5.515 


17 


129940 


1913501 


2668241 


4338715 


12.500 


18 


279415 




7652212 


12785056 


27.109 


19 


601742 




22013683 


37739184 


72.579 


20 


1297671 




63497798 


111562926 


163.422 


21 


2802318 




183589726 


330215133 


383.640 



Table 1. Enumeration of conjugacy classes to generic level = 21, 
with run times measured in seconds. The tree numbers were computed 
in [DM]. 



Levels 17/16 


Levels 18/17 


Levels 19 / 18 


Levels 20 / 19 


Levels 21/20 


2.941 


2.947 


2.952 


2.956 


2.960 



Table 2. Computing the growth: ratios of numbers of conjugacy 
classes at consecutive levels. 



classes associated to each truncated spine. The ideas and proofs follow the treatment 
of cubic polynomials in [BHl] . jBH2] . [Br], and jBDK] . 

1.1. Results of the computation. An implementation of the algorithm was written 
with Java. We compiled the output in Table[T]to level = 21, together with run times 
(Processor: 2.39GHz Intel Core 2 Duo, Memory: 2 x 1GB PC2100 DDR 266MHz). 

1.2. The tree of cubic polynomials. We can define a tree Tconj of conjugacy 
classes of cubic polynomials in the shift locus as follows. For each generic level 
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Figure 1. The tree Tconj of conjugacy classes to generic level = 5. 

> 1, let V^(A^) be a set of vertices consisting of one vertex for each topological 
conjugacy class. We connect a vertex v in V{N + 1) to a vertex w in V^(A^) by 
an edge if for each representive f of v and g of w, there exists an e > so that 
the restrictions f\{Gf > (G/(ci)/3^-i) - e} and g\{Gg > {Gg{ci)/3^-^) - e} are 
topologically conjugate. In Figure [l| we have drawn this tree %onj to level = 5. 
The growth rate of the number of conjugacy classes as iV — > cx) corresponds to a 
computation of the "entropy" of this tree. 

The tree %onj can also be constructed in the following way. Let MP3 denote 
the space of conformal conjugacy classes of cubic polynomials; it is a two-dimensional 
complex orbifold with underlying manifold isomorphic to . Every cubic polynomial 
is conformally conjugate to one of the form 

f{z) = + az + b, 

which can be represented in MP3 by (a, 6^) G C^. 
The critical escape-rate map 

g : MP3 ^ 

is defined by Q{f) = {Gf{ci), Gf{c2)) where the critical points {ci, C2} of / are labeled 
so that Gf{ci) > G/(c2); it is continuous and proper |BHlj . The fiber of Q over the 
origin in is the connectedness locus C3, the set of polynomials with connected Julia 
set. If we restrict Q to its complement MP3 \C3, there is an induced projectivization: 

^:MP3\C3^[0,1] 

defined by / Gf{c2)/Gf{ci). The quotient space of MP3 \C3 formed by collapsing 
connected components of fibers of Q to points is a (completed) tree P7^*: over (0, 1] 
it forms a locally finite simplicial tree while ^~^(0) forms its space of ends. 

In |DP2j it is proved that the edges of the tree P7^* correspond to generic topological 
conjugacy classes. Thus, the combinatorial tree Tconj is dual to the shift locus tree in 
PT^*; each edge in P7^* corresponds to a vertex in Tconj- 
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The tree FT^ comes equipped with a projection to the space of cubic trees F% 
introduced in |DMj . In |DMj the growth of the number of edges in FT^ was stud- 
ied (see Table [l| third column), but the value of the entropy was left as an open 
question. Furthermore, the tree Tconj is a quotient of the tree of marker autorphisms 
A^3 introduced in |BDK] : the quotient is by the monodromy action from a twisting 
deformation (see |Brj ). The entropy of A^a was easily shown to be log 3, and so the 
entropy of T^onj (or equivalently, of P7^*) is no more than log 3. 

Question. Let (p{N) denote the number of vertices in Tconj at level N. Is 

lim — \og(f{N) = log 3 ? 

In Table |2| we show the ratios of the number of conjugacy classes in consecutive 
levels. As the level increases, the computed ratios increase, conjecturally limiting on 
3. 



2. The r functions 

2.1. The r-function of a polynomial. Fix a cubic polynomial / with disconnected 
Juha set, and let ci and C2 be its critical points, labeled so that Gf (02) < G'/(ci). 
For each integer n such that Gf{c2) < Gf{ci)/3'^~^, we define the critical puzzle piece 
Pn{f) as the connected component of {z : Gf{z) < G'/(ci)/3"~^} containing C2, and 
set 

Tfin) = max{j < n : /""^(ca) G P,(/)}. 

Recall that the tableau or marked grid of f is an array {Mf{j, k) G {0, 1} : j,k > 0}, 
defined by the condition 

Mf{j,k) = i ^ f(^c2)eP,U-)- 

We depict a marked grid as a subset of the 4th quadrant of the Z^-lattice, where j > 
represents the distance along the negative y-axis and A; > respresents the distance 
along the positive x-axis. The values of Tf can be read directly from the marked grid: 
beginning with Mj(n, 0) = 1, Tf{n) is the j-coordinate at the first non-zero entry 
when proceeding "northeast" from (n, 0). In fact, the orbit {Tj(n) : A; > 0} consists 
of the j-coordinates of all non-zero entries along the diagonal Mf{n — i,i). Thus, the 
marked grid can be recovered from the r-function by: 

r 1 if J = A; = 

Mf{j, k) = < 1 if j = T™(j + k) for some m > 
otherwise 

Branner and Hubbard |BH2l Theorem 4.1] showed that marked grids associated to 
cubic polynomials are characterized by a simple set of rules. A marked grid of size 
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(which may be infinite) is an array {M{j, k) G {0, 1} : j, A; > and j + k < N} 
which satisfies the following rules: 

(MO) For each n<N, M(n, 0) = M(0,n) = 1. 

(Ml) If M(j, k) = 1, then M{1, k) = 1 for all I < j. 

(M2) If M(j, k) = 1, then M{j -i,k + i) = M{j - i, i) for all < i < j. 

(M3) If j + A; < A^, M(j, k) = 1, M(j + 1, A;) = 0, M{j - = hr < i < m, 

and M{j — m + 1, m) = 1, then M{j — m + l,k + m) = 0. 
(M4) Ifj + k < N, M{j,k) = 1, M(l,j) = 0, M(j + 1,A;) = 1, and M{j -i,k + i) = 

for all < i < j, then M(l, j + k) =0. 

The rule (M4) was omitted in |BH2] , though it is necessary for their proof. It appears 
as stated here in [Kil Proposition 4.5]; an equivalent formulation (in the language of 
r- functions) was given in |DM] . 



2.2. Properties of tau-functions. Let N denote the positive integers {1, 2, 3, . . .}. 

We consider the following five properties of functions r : N ^ N U {0}. 

(A) r(l)=0 

(B) T{n + 1) < T{n) + 1 

From (A) and (B), it follows that T{n) < n for all n G N; consequently, there exists 
a unique integer ord(n) such that the iterate T"'^'^^^\n) = 0. 

(C) If r(n + 1) < r^(n) + 1 for some < k < ord(n), then r(n + 1) < r'''+^(n) + 1. 

(D) If r(n + 1) < T''{n) + 1 for some < k < ord(n), and if T{T^{n) + 1) = 
r^+i(n) + 1, then r(n + 1) < r^'+^(n) + 1. 

(E) If ord(n) > 1 and oid{T°"^^''^-^ {n) + 1) = 1, then r(n + 1) ^ 0. 

Proposition 2.1. For any positive integer N , a function 

r:{l,2,3,...,Ar}^NU{0} 

or a function 

r : N^NU{0} 

is the T-function of a cubic polynomial if and only if it satisfies properties (A)-(E). 

We say the r-function is admissible if it satisfies properties (A)-(E). 

The proof is by induction on A^. It is not hard to see that the r function must satisfy 
these rules, by doing a translation of the tableau rules. Conversely, any tau function 
satisfying properties (A)-(E) determines a marked grid satisfying the 4 tableau rules. 
Property (E) is another formulation of the "missing tableau rule" (M4) appearing in 
[Ki] and [DM]. 
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2.3. Algorithm to inductively produce all r-functions. If a r-function has do- 
main {1, 2, 3, ... , N}, we say it has length N. The markers of a r with length are 
the integers 

{m e {1, . . . , N - 1} : T{m + 1) < r(m) + 1}. 
Let k be the number of markers which appear in the orbit 

N ^ r{N) ^ . . . ^ r^'^^i^XN) = 0, 
and label these k markers by l[, /g, . . . , so that 

N = 1'q>1[>1'^> ■■■>l[>0. 
For each < i < k, let k = t^Q so that 

t{N) = lo>k> ■■■>lk>0. 

Theorem 2.2. Given an admissible r-function of length N , an extension to length 
N + 1 is admissible if and only if 

t{N + 1) = /j + 1 for some <i<k 

or t{N +l)=0iflk>0ork = 0. 



Proof. The theorem follows from Proposition 2.1 Property (C) implies that t{N + 1) 
must be either or of the form t^{N) + 1. Property (D) implies that t{N + 1) must 
be either or of the form h + 1. Property (E) implies that t{N + 1) 7^0iffc>0 
and Ik = 0. On the other hand, if = 0, both t{N + 1) = /q + 1 = r{N) + 1 and 
r(A^ + 1) = are admissible. □ 



3. The truncated spine 

Suppose / is a cubic polynomial of generic level N. Introduced in |DPlj . the 
truncated spine of / is a combinatorial object which carries more information than 
the r-sequence though it does not generally determine the topological conjugacy class. 
(It determines the tree of local models for /, studied in |DP1] .) Here, we describe how 
to inductively construct truncated spines, and we compute the number of extensions 
to a truncated spine of length + 1 from one of length A^. We show exactly how 
many distinct extensions correspond to a choice of r-function extension. 

3.1. The truncated spine of a polynomial. Fix a cubic polynomial / of generic 
level A^. The truncated spine is a sequence of A^ finite hyperbolic laminations, one 
for each connected component of the critical level sets of Gf separating the critical 
value /(ci) from the critical point C2, together with a labeling by integers < A^. 

Specifically, beginning with the level G/(ci) of fastest-escaping critical point, we 
identify the level set {Gf = Gf{ci)} with the quotient of a metrized circle. The level 
curve is topologically a figure 8, metrized by its external angles, giving it total length 
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27r. The critical point Ci lies at the singular point of the figure 8, identifying points 
at distance 27r/3 along a metrized circle. Thus the associated hyperbolic lamination 
in the unit disk consists of a single hyperbolic geodesic joining two boundary points 
at distance 27r/3. The lamination is only determined up to rotation. We mark the 
complementary component in the disk with boundary length 47r/3 with a to indicate 
the component of {Gf < G/(ci)} containing the second critical point C2- 

For each critical level set {Gf = G/(ci)/3"}, < n < A^, we consider only the 
connected component which separates C2 from ci. External angles can be used to 
metrize the curve, normalizing the angles so the connected curve has total length 27r 
for each n. The curve can thus be represented as a quotient of a metrized circle; the 
associated hyperbolic lamination consists of a hyperbolic geodesic joining each pair 
of identified points. The gaps of the lamination are the connected components of 
the complement of the lamination in the disk; each gap corresponds to a connected 
component of {Gf < G/(ci)/3"}. A gap in a lamination is labeled by the integer 
A; > if the corresponding component contains f^{c2)- 

In Figure [2| we provide examples of truncated spines for two polynomials of generic 
level N = 5 with the same r-sequence. 

Two truncated spines of length N are equivalent if for each < n < A^, the labeled 
laminations at level n are the same, up to a rotation. In particular, the labeling of the 
gaps must coincide. It is shown in |DPlj that a truncated spine (together with the 
heights {Gf{ci), Gf{c2)) of the critical points) determines the full tree of local models 
for a polynomial of generic level A^. In particular, it carries more information than 
the tree of /, introduced in |DMj . and its r-function; it does not, however, determine 
the topological conjugacy class of the polynomial. 

The central component of a lamination is the gap containing the symbol 0; it cor- 
responds to the component containing the critical point C2. All non-central gaps are 
called side components. 

3.2. Prom truncated spine to r-function. Fix a truncated spine of length A^. 
Recall that every r-function satisfies r(l) = 0. For A^ > 1 and each < n < A^, we 
can read r(n) directly from the truncated spine by: 



To compute t{N), we consider the set 

L{N) = {j : the central component in lamination j is labeled by (A^ — j — 1)}- 
We then have 



r(ra) = max{j : lamination j is labeled by {n — j)}. 




1 + max{j : j E L{N)} ii L{N) ^ $ 
if L{N) = 
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Figure 2. Two truncated spines associated to polynomials of generic 
level = 5 with the same r-sequence 0, 1, 0, 1, 0. 

3.3. Extending a spine of length A^. Fix a truncated spine of length N. It follows 
directly from the definitions that an extension to length + 1 is completely deter- 
mined by the location of the label N — t{N) at level t{N). Any choice of central or 
side component is admissible: it determines the local model for an extended tree of 
local models (from |DP1] ). 

The lamination at level A^ is then constructed by taking a degree 2 branched cover 
of the lamination at level t{N) branched over the gap containing the label N — t{N). 
See |DPlj for a general treatment of branched covers of laminations. The labels are 
added inductively: for each iterate 1 < n < ord(A^), the label (A^ — r"(A^)) is placed 
in the gap at level t^{N) which is the image of the gap containing r"^^(A^) at level 
r"-i(Ar). 

3.4. Computing the number of extended spines for each choice of t{N + 1). 



Fix a truncated spine with its r-function of length A^. As in ^:2.3, the markers of r 
are the integers 

{m e {1, . . . , N - 1} : T{m + 1) < r(m) + 1}. 
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The marked levels of r are all integers in the forward orbits of the markers: 

{/ > : / = T'^{m) for marker m and n > 0} U {0}; 

we say is marked even if there are no markers. 

As before, we let k be the number of markers which appear in the orbit 

N ^ t{N) ^ ...^ r°"^(^\N) = 0. 

Label these k markers hj l[,l2y ■ ■ ■ y^'k that 

N = l'o>l[>l2> ■■■>l'k>0. 

For each < i < k, let k = t{1[) so that 

r(iV) = k>h> ■■■>lk>Q. 

For each < i < fc, define be the condition that 

and define Uk so that T"''^{lk) = 0. For < i < j < /c + 1, we set 



1 if r(/^ + 1) = + 1 
otherwise 



where by convention we take Ik+i = — 1. Note that r(/^ + 1) = for every r, so 
6{k,k + l) = 1. 

At level Iq = t{N), there are 2°^^^'-''^ = 2"o+"i+'''+"'= side components and one central 
component. Labeling the central component with the integer (A^ — t{N)) uniquely 
corresponds to the choice of t{N + 1) = t{N) + 1. For each i > 0, the number of side 
components which correspond to the choice t{N + 1) = /j + 1 is 

2no(2"i(2-2(...(2"-i )~5{2,i))-5{l,i)); 

as above, we take Ik+i = —1. It remains to consider how many distinct truncated 
spines these side components determine. 
The symmetry of r is 

s = min{n > : t"(/o) is a marked level}. 



Note that s < uq. To each admissible choice for t{N + 1) (from Theorem 2.2) we 
define the (A^ + l)-th spine factor of r. If t{N + 1) = /j + 1, then 

/ 2"o-(2"^(2"^(---(2"'--5(z-l,0)----)-5(2,z))-5(l,z)) if2>0 
SF(r,Ar+l) = | ^ .^^^^ 

where, as above, we take Ik+i = — 1. It is now straightforward to see: 
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Theorem 3.1. For any t -function of length N, the number of truncated spines with 
this T-function is: 

N 

Spines (r) = JJ SF (r,j). 



4. Complete algorithm 

We combine the resuhs of the previous sections to produce an algorithm for the 
complete count of topological conjugacy classes for cubic polynomials with generic 
level A^. 

Fix a r-function of length N. Recall that the markers of r are the integers 
{me{l,...,N -1} : r(m + 1) < r(m) + 1}, 
and the marked levels are: 

{/ > : / = T^{m) for some n > and marker m} U {0}; 
Let L be the number of non-zero marked levels. 



For each n < N, the order of n was defined in ^2.2; it satisfies r°^^^^\n) = 0. For 
each marked level /, compute 

I 

mod(0 = J]2-"'^« 

and 

t{l) = min{n > : nmod(/) G N}. 

The quantity mod(/) represents the sum of relative moduli of annuli down to level /, 
while t{l) is the numer of twists required to return that marked level to its original 
configuration. We define 

T(r) = max{t(/) : / is a marked level} 

or set T(r) = 1 if r has no marked levels. The twist factor is defined by 

2^ 

TFfr) 



T(r) 



From Theorem 3J^, the number of truncated spines with this r-function is: 

N 

Spines(r) = JJSF(r, j). 

By |DPlj . the number of topological conjugacy classes associated to r is then 

Top(r) = Spines(r) ■ TF(r). 
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Combining these computations with Theorem 2.2, it is straighforward to automate an 
inductive construction of all r-functions of length A^, and we obtain an enumeration 
of all truncated spines and all topological conjugacy classes of generic level N. 
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